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Method and device for using & multi-channel measurement signal 
in determining the current distribution of an object 

FIELD OF THE INVENTION 

The present invention relates to a new and 
advanced method for using a multi-channel measurement 
signal in source modelling. In particular , the present 
invention relates to a new way of converting the meas- 
urement signals measured using a multi-channel measur- 
ing device that measures an irrotational and source- 
less vector field into a form optimal from the stand- 
point of the source modelling. 

BACKGROUND OF THE INVENTION 

The electrical operation of an object to be 
examined can be examined e.g. by measuring the mag- 
netic field produced by the electric currents of the 
object using sensors placed outside the object. The 
modelling of a source distribution thus obtained based 
on the measured signals is, however, very difficult 
because each magnetic field distribution can be pro- 
duced by many different source distributions. In other 
words, a source distribution cannot be solved unambi- 
guously based on the measured signals, so to solve the 
problem, different restrictive conditions must be set, 
such as some parametric model based on prior informa- 
tion for a current, or a non-parametrised norm re- 
striction. 

For the non-parametrised modelling of a con- 
tinuous current distribution, minimum norm estimates 
are usually used, in which there is an attempt to ex- 
plain a signal measured using a multi-channel measure- 
ment device with a current distribution whose norm is 
as small as possible. As the norm, usually an LI or L2 
norm is selected, of which the previous is the sum of 
the lengths of the current elements over the selected 
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volume, and the latter is the sum of the squares of 
the lengths of the current elements over the selected 
volume. The calculation of the minimum norm estimates 
has been described e.g. in publications ^Interpreting 
magnetic fields of the brain: minimum norm estimates", 
M.S. Hamaiainen et al, Medical & Biological Engineer- 
ing & Computing, Vol. 32, pp. 35-42, 1994, as well as 
"Visualization of magnetoencephalographic data using 
minimum current estimates", Uutela K. et al, Neurol- 
mage, Vol. 10, pp. 173-180, 1999. 

Conventional minimum norm estimates involve 
inherent problems such as slowness of calculation and 
susceptibility to noise. For example, in the case of 
an L2 norm, one needs an inverse matrix of matrix G, 
whose element (i, j) contains the inner product of the 
switching fields of the ith and jth measurement sen- 
sor, so one must calculate these inner products for 
each pair of sensors. The switching field is so deter- 
mined that the signal measured by a sensor is the pro- 
jection of the current distribution for the switching 
field of the sensor in question. The noise problems 
are due to the fact that matrix G calculated for the 
sensors is susceptible to noise, so in the calculation 
of its inverse matrix, regularisation is needed in. the 
practical situations. 

Regularisation methods, such as the break-off 
regularisation of the singular value decomposition, 
usually are non-intuitive, and usually also to be 
solved for each case specifically. A regularisation of 
the wrong type may lead to an erroneous modelling re- 
sult. 

Therefore, source modelling nowadays still 
involve problems such as the hardness and slowness of 
the computation, the possible errors caused by noise, 
as well as the case-specificity due to the regularisa- 
tion. Further, as stated above, the regularisation may 
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cause considerably errors to the final computation re- 
sult. 

OBJECTIVE OP THE INVENTION 

The objective of the invention is to elimi- 
nate the drawbacks referred to above, or at least sig- 
nificantly to alleviate them. One specific objective 
of the invention is to disclose a new type of method 
which can be used to considerably lighten and acceler- 
ate the computation associated with the modelling of a 
continuous current distribution, as well as to lessen 
the problems with noise. 

As for the features characteristic of the in- 
vention, reference is made to them in the claims. 

DESCRIPTION OF THE INVENTION 

The present invention relates to a new kind 
of manner of determining the continuous current dis- 
tribution of an object being examined using basis vec- 
tor components of signal space calculated from the 
measured signals. The components in question have been 
so selected that they describe the features, as inde- 
pendent as possible, of the current distribution being 
examined, which enhances the computation and makes it 
more accurate. 

The basic idea of the invention is that be- 
cause the computation of the inner products of the 
sensor fields is hard and difficult using a conven- 
tional set of sensors, it is worth using a special set 
of sensors, whose switching fields are orthogonal and, 
if possible, to be analytically computed. In princi- 
ple, this can be implemented by a suitable physical 
set of sensors. As a suitable physical set of sensors 
is, however, often quite difficult to manufacture, it 
is, in most of the cases, advantageous to use virtual 
sensors computationally generated from a conventional 
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set of sensors, i.e. the measurement signals are. con- 
verted into other ones by a suitable conversion so 
that they correspond to the signals that the virtual 
measurement device would have measured. At the same 
time, it is possible, if necessary, to eliminate the 
signals associated with external interference. This 
conversion has been described e.g. in patent applica- 
tion FI2003Q392, which is incorporated herein by ref- 
erence. After the conversion, the source modelling is 
performed in an optimal manner using the basis vector 
components of the signal space instead of the actual 
measurement signals. One substantial feature of the 
invention is that after the conversion, the source 
model need not be any more regularised. 

Thus, the present invention relates to a 
method for determining the current distribution of an 
object by measuring the magnetic fields in the vicin- 
ity of the object using a multi-channel measurement 
device. Advantageously, at least one measurement sen- 
sor corresponds to each channel, and the object is ap- 
proximated by a spherical-symmetrical conductor. The 
object can be e.g. a human being's head. 

According to the invention, a multi-channel 
measurement signal corresponding to each measurement 
sensor is converted into the signals of a predeter- 
mined set of virtual sensors, and the current distri- 
bution of the object being examined is determined by 
depth r from the signals of the set of virtual sensors 
in a beforehand selected orthonormal function basis. 
In that case, the estimation of a current distribution 
is fast and robust. Further, to achieve the set of 
signals corresponding to the set of virtual sensors, a 
multi-pole development is calculated from a multi- 
channel measurement signal. A multi-pole development 
can be calculated in two ways: by taking into account 
the magnetic fields outside the object being measured, 
or by ignoring them. 
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Advantageously, as the orthonormal function 
basis, a basis with the following form is selected: 

M = lta(r)im(M. 

wherein f ( r) is a radial function to be 
freely selected and Jc^i 0 *?) is so-called vector 
spherical harmonic. In that case, it is . possible to 
place the function basis into a current distribution 
equation, and the coefficients of the current distri- 
bution are analytically solved based on the equation: 

wherein y t is a constant associated with or- 
der 1 and R is the radius of the sphere to be exam- 
ined. Advantageously, function f(r) is used to adjust 
the depth weighing of a current distribution model. 

Furthermore, the invention relates to a meas- 
urement device for determining the current distribu- 
tion by measuring magnetic fields in the vicinity of 
the object* The measurement device includes a set of 
measurement channels that measure an irrotational and 
sourceless vector field, whereby at least one measure- 
ment sensor corresponds to each channel, and process- 
ing means for processing the measurement signal. Ad- 
vantageously, the object is approximated with a 
spherical-symmetrical conductor. 

According to the invention, the processing 
means include a conversion module for converting a 
multi-channel measurement signal corresponding to each 
measurement sensor into the signal of a predetermined 
set of virtual sensors; and calculation means for de- 
termining the current distribution of an object being 
examined or for calculating by depth r from the sig- 



WO 2005/030051 



PC17FI2004/000532 



6 

nals of the set of virtual sensors in a beforehand se- 
lected orthonormal function basis. In one embodiment, 
the calculation means are arranged to calculate a 
multi-pole development from a multi-channel measure- 
ment signal. 

The invention enables one to considerably 
lighten and accelerate the calculation associated with 
the modelling of a continuous current distribution. 
The invention further enables one to reduce the prob- 
lems associated with noise. Further , the invention 
simplifies the regularisation of a source model, or 
eliminates the need for it, thus also significantly 
reducing the possibility of error. 

LIST OF FIGURES 

In the following, the invention will be de- 
scribed in detail by means of examples with reference 
to the accompanying drawing, in which 

Fig. 1 represents one measurement device in 
accordance with the present invention; and 

Fig. 2 is a flow chart illustrating one 
method in accordance with the present invention. 

DETAILED DESCRIPTION OF THE INVENTION 

Fig. 1 represents one measurement device in 
accordance with the present invention comprising a set 
of measurement channels l 1 , l 2 ,...l 5 that measure an 
irrotational and sourceless vector field, in which 
case at least one measurement sensor 2 1 , 2 2 , ...2 5 cor- 
responds to each channel; and processing means 3 for 
processing the measurement signal. Advantageously, the 
processing means have been implemented by means of a 
computer. Further, the processing means include a con- 
version module 4 for converting a multi-channel meas- 
urement signal corresponding to each measurement sen- 
sor into the signals of a predetermined set of virtual 
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sensors, and calculation means 5 for determining the 
current distribution of an object being examined. 

Further Fig. 1 shows the object being meas- 
ured K, into the vicinity of which the measurement 
sensors 2 1 , 2 2 , ...2 s have been placed. Inside the ob- 
ject being measured there is presented a current loop 
which describes the source of the magnetic field being 
measured. The object being measured can be e.g. a hu- 
man head, and there can be several current sources. 

Fig. 2 represents the main steps of one pre- 
ferred embodiment of the present invention. At first, 
a set of sensors is used to measure a multi-channel 
measurement signal, step 21. After that, the signal is 
converted to correspond to the signal measured using a 
so-called set of virtual sensors, step 22, making the 
mathematical operations simpler. Finally, it is possi- 
ble to simply calculate the current distribution in an 
object from the signal of the set of virtual sensors, 
step 23, i.e. in practice, to describe the places and 
strengths of the current loops inside a sphere or 
head. 

In the following section, the mathematical 
background and grounds of the invention are described. 
When the magnetic fields are converted into coeffi- 
cients M &l =a te +/fe to associated with the basic solution 
r -v+»Yim(P,<p) of the Laplace equation, wherein i is an 
imaginary unit, they can be expressed by means of the 
current distribution in spherical coordinates 

(r 9 0 9 <p), whereby they are of the form: 

tot* =ri\r l $to{0,<P)*$(P)dv (1) , 

wherein the integration is performed over the 
volume being examined, y t is a constant associated 

with order 1 and Jhm(0 9 <p) is so-called vector spherical 
harmonic. 
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This form is apparent e.g. from publication 
"Multipole expansions of electromagnetic fields using 
Debye potentials" , C.G. Gray, American Journal of 
Physics, Vol. 46, pp. 169-179, 1978. The expression 
mentioned above is of the switching field form, 
wherein the switching field of the multi-pole coeffi- 
cient Mim is of the form: 

4(ft-r^ ta (^) (2). 

On the other hand, the vector spherical har- 
monics form by depth r an orthonormal basis, so with 
the depth in question, the current distribution can be 
presented in the function basis in question: 

M = 2 £c^(r)&*(0 9 p) ( 3 ) , 

wherein f (r) is some radical function. 

When as the volume to be examined, spherical 
volume is selected, by placing the previous, expression 
into the equation (1), the coefficients of the current 
distribution can be analytically solved: 

c^mj^r^mdr^ (4), 

wherein fx is a constant associated with or- 
der 1 and J? is the radius of the sphere to be exam- 
ined. The previous equation (4) indicates that the co- 
efficients of a current distribution model presented 
in an orthonormal basis can be solved based on coeffi- 
cients i^Lm in a completely trivial manner using ana- 
lytical expressions without any kind of regularisa- 
tion. This is computationally very fast and numeri- 
cally stabile. Function f(r) is freely selectable, and 
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can be used to adjust the depth weighing of a current 
distribution model. 

Furthermore, it must be noted that in the 
case of a spherical conductor, it is possible to use 
as basis functions also some other orthogonal basis or 
a basis whose inner products can be otherwise quickly 
calculated. This is achieved e.g. by slightly breaking 
the orthogonality in a manner known per se, or by us- 
ing a basis which is not orthogonal but whose inner 
products can be calculated beforehand. 

The invention is not limited merely to em- 
bodiment examples referred to above, instead many 
variations are possible within the scope of the inven- 
tive idea defined by the claims. 



